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Abstract 

GNSS-R signals contain a coherent and an incoherent component. A new algorithm for 
coherent phase altimetry over rough ocean surfaces, called PARFAIT, has been developed and 
implemented in Starlab's STARLIGHT 1 GNSS-R software package. In this paper we report 
our extraction and analysis of the coherent component of LI GPS-R signals collected during 
the ESTEC Bridge 2 experimental campaign using this technique. The altimetric results have 
been compared with a GPS-buoy calibrated tide model with a resulting precision of the order 
1 cm. 

Keywords: Passive radar, GNSS, GPS, Galileo, GNSS-R, GPS-R, altimetry, PARIS, PIP, 
PARFAIT, coastal applications. 



1. Introduction 

Specular reflections dominate medium to 
short wavelength electromagnetic forward 
scattering on the ocean, examples of which in- 
clude GNSS and solar reflections. As reported 
in during the last decade many GPS- 

R (Global Positioning System Reflections) ex- 
perimental campaigns have now been success- 
fully carried out. A partial and surely incom- 
plete list is provided in Table 1. In this paper 
we focus on the potential of GNSS-R (Global 
Navigation Satellite System Reflections) for al- 
timetric coastal applications. The techniques 
developed, however, can also be implemented 
in other scenarios (airborne, spaceborne). 

The specularly scattered field is composed 
of a coherent component and a random, Hoyt- 
distributed component ;2]. When the surface 
is very rough, the latter becomes incoherent 
and the former becomes very small. In fact, if 
the surface height distribution is normal with 
deviation o^, then 

(r 2 ) ^ n 2 e -(4™ c co S e/A) 2 + (1) 

n(l-e-( 4 ^- cos9 / A > 2 ), 

where (r 2 ) is the power average, n is the num- 
ber of scatterers, A the EM wavelenght and 
the local incidence angle \2'A\ I24| . 

GNSS-R signals thus contain a coherent 
and an incoherent component. In companion 
papers we discuss the analysis of the incoher- 
ent component for sea state monitoring |15II25| 
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and for code altimetry [37J . Here we present a 
new approach for the extraction and analysis 
of the coherent component of GNSS reflected 
signals to perform phase altimetry. 

The data discussed here was collected by 
ESA/ESTEC during the Bridge-2 experiment. 
The experiment aimed at gathering direct and 
reflected GPS signals from antennas located 
about 18 m above the mean sea level of an 
Estuary in the North sea of Holland. For 
more information on the experimental setup, 
the reader is directed to |21| . 

This paper is structured as follows: 

• Discussion on coherence properties of re- 
flected signals and their use for phase al- 
timetry. 

• Analysis of the direct and reflected sig- 
nals and EM field extraction. 

• Implementation of PARFAIT altimetry. 

• Comparison with other data and discus- 
sion of the altimetric results. 

2. The complex field 

The importance of retrieving the coher- 
ent part of the EM field backscattered by the 
sea surface stems from its altimetric content. 
Measuring the phase of the coherent compo- 
nent allows for accurately estimating the delay 
of the reflected signal with respect to the the 
direct one, i.e., for estimation of the temporal 
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Main 
Author 


Inst. 


Date 


Notes 


Garrison 


El 


JPL 


1996 


A normal GPS receiver was used. Demon- 
strated reception/tracking of reflected signals 
over relatively calm waters. Concluded that 
more complex receiver is needed. 


Martin-Neira 
Caparrini 


[HE] 


ESA 


1997 


GNSS-R PARIS altimetric experiment from a 
Zeeland bridge. C/A code used for correla- 
tion leading to an altimetric accuracy in the 
order of 3 meters after 1 second (1% of the 
chip length). 


Garrison 




JPL 


1997 


Widening of correlation function in rough seas 
demonstrated. Application for sea state from 
air. 


Komjathy 


16 


CCAR 


1998 


Aircraft experiments, 3-5 km altitude. 


LaBrecque 


EE 


NASA 


1998 


The first spaceborne observation of GPS sig- 
nals reflected from the ocean surface. 


Cardellach 

Ruffini 

Garrison 




IEEC 


1999 


Balloon experiment. Successful detection of 
reflected signals at 38 km of height with low 
gain antenna. Sea surface winds retrieved 
with ~2 m/s error. 


Cardellach 
Ruffini 


n 


IEEC 


1999 


First ESA aircraft experiment. Some sig- 
nals detected, DDM produced, but experi- 
ment failed due to hardware problems. 


Armatys 


in 


CCAR 


2000 


Wind speed and directions obtained from re- 
flected GPS signals are compared to the Sea- 
Winds scatterometer on-board QuikSCAT. 


Garrison 




JPL 


2000 


With GPS-R airborne data, retrieval of the 
wind speed with a bias of less than 0.1 m/s 
and with a standard deviation of 1.3 m/s. 


Zavorotny 


m 


CCAR 


2000 


Fundamental theoretical work. Comparison of 
experimental and theoretical waveforms. 


Zuffada 


m 


JPL 


2000 


Lakeside experiment, with an almost flat sur- 
face (no roughness). Centimetric phase al- 
timetry. 


Martin-Neira 
Ruffini, 
Serra, 
Col- 

menarejo 


m 


ESA 


2000 


The pond experiment was designed to test 
some key issues in the PARIS Interferomet- 
ric Processor (PIP) concept. The PIP con- 
cept is based on the use of dual-frequency car- 
rier measurements to exploit the correlations 
in the scattered signals at similar frequencies. 


Ruffini 
Caparrini 


□ 


Star- 
lab 
IEEC 


2001 


GPS-R LI data collected from the Casablanca 
drilling platform by IEEC has been analysed 
at Starlab. 


Martin-Neira 


m 


ESA 


2001 


The experimental campaign which is the ob- 
ject of this work. 


Cardellach 

Starlab 

Team 




IEEC/ 
Star- 
lab 


2001 


GPS-R data collection from airborne plat- 
form. Campaign performed within the 
ESA/ESTEC project "PARIS Alpha". Data 
processed under ESA/ESTEC projects 
"PARIS Alpha" and "OPPSCAT 2". 


Starlab 
Team 


[E3[27j 


Starlab 


2002 


GPS-R data collection from airborne plat- 
form. Retrieval of altimetric profile matching 
Jason- 1. (ESA/ESTEC "PARIS Gamma") 


Starlab 
Team 


m 


Starlab 


2003 


GPS-R data collection from Barcelona Har- 
bour (HOPE campaign, Starlab Oceanpal 
project). 
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lapse. This is the essential measurement for 
altimetry |26| . 

In order to collect the complex EM field, 
the complex signal is generated from the 
real one. Although this operation is often 
performed by the receiver front-end, in the 
Bridge-2 experiment only the in-phase com- 
ponent of the signal was sampled at high 
frequency and stored on digital tape. The 
quadrature component was generated after- 
wards. The process is illustrated in Figure (QJ. 
We can then represent the direct signal re- 
ceived at the antenna input 2 as 

S&) = A d ■ C(t) ■ D(t) ■ e *(""+^)' + Vd , (2) 

where Ad is the direct signal amplitude, C(i) 
represents the C/A code, D(t) the naviga- 
tion code, ojli the LI carrier frequency, tod 
the Doppler frequency offset, and rjd (thermal) 
noise. The reflected signal at low altitudes can 
be modelled by 

S r (t) = C(t) ■ D{t) • e '(^i+^)t. ( 3 ) 

(A, ■ e 2mC ' X + 0(tj) + V r, 

where A r is the reflected signal mean ampli- 
tude, 0(t) represents the perturbation due to 
ocean motion and C the reflected signal extra 
path length. In coastal applications 0(t) is 
a relatively slowly varying quantity with zero 
mean, while C, which contains the geophysical 
tide signal, can be considered effectively frozen 
during correlation processing. 

After modulation with a local oscillator of 
frequency u>li—ojif and low-pass filtering, the 
signal will have a residual carrier at uid+uiF- 
This signal is mixed with a phasor at fre- 
quency ujjf + Qd, where u>d is an estimate of 
the Doppler frequency for the satellite under 
investigation, and finally low-pass filtered. 

With the assumption that the navigation 
bit is constant during the integration time 
(which is correct if the correlation is bit- 
aligned and the coherent integration time Tg 
is less than 20 ms), and considering that dur- 
ing an integration time interval the value of 
Aoj d is constant, the complex p-th sample of 
the correlation coefficient for the direct signal 

2 This contains only the C/A code part. The P code 
correlation of the signal with replicas of the C/A code- 



writes 

C p ~ l -A d D k R p e- lA ^ T ^-^ ■ 



T, 0111 I 2 & I 

sm [s^T' ) 

where T s is the sampling interval and R p the 
corresponding correlation coefficient function. 
For the reflected signal we can write an equiva- 
lent expression, modulated by the slowly vary- 
ing phasor A r -exp(2iri£/A)+0(t). For coastal 
applications we can assume there will be little 
filtering of this quantity by the coherent in- 
tegration process, as the ocean moves slowly 
compared to coherent integration times (a few 
ms). 

In the case of the direct signal, we can eas- 
ily track the carrier phase. To this end, the 
delta-phases obtainable from equation J3J can 
be accumulated using 

P +i -4>v = 

= Im(k>g^|i) = -Au dp T E . (5) 

This equation holds while Au>d p+1 ~ Aujd p - 
This is a good approximation, since the time 
during which this variation is measured is the 
coherent integration time. 

The main advantage of using this algo- 
rithm for phase tracking is that, due to its dif- 
ferencing nature, it allows for easy detection 
of the navigation bit 7r radians phase change. 

Figures [3 to |H| which illustrate these con- 
cepts, refer to the processing of another set of 
GPS-R data — collected during the Casablanca 
oil platform Experiment. This Repsol owned 
drilling platform is about 40 km off the coast 
of Tarragona, Spain (40°43'4"iV, 1°21'34"£). 
The measurement campaign took place on 
March 14th, 2000. 

In Figure the histogram of the delta- 
phases is shown. The x-axis represents cycles 
and the y-axis is in arbitrary units. Most 6- 
phase values are clearly concentrated around 
zero. Other values gather just before ±tt. 
These values represent in fact small values to 
which ±7r radians have been added on occur- 
rence of a navigation bit transition. 

In Figure El the direct signal phase with 
and without navigation bit correction is plot- 
ted. In Figure 01 the phase for the (navi- 

component can be neglected thanks to the subsequent 
-the two codes are orthogonal. 
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gation bit corrected) direct and reflected sig- 
nals is shown. The effect of the reflection on 
the sea surface is clearly visible in the large 
variations present in Figure |4(b)| with respect 
to Figure 4(a) In Figure |51 and El the am- 
plitude magnitude and the complex vector of 
the direct and reflected fields, respectively, are 
shown. In Figure |H1 a simulation of the LI 
GPS complex field phasor dustball after reflec- 
tion, akin to the one in Figure 6(b) is shown. 



The simulation parameters have been chosen 
to match the Casablanca experiment sea state. 
Those were reported as a "quite calm sea with 
a gentle breeze", with SWH of about 0.7 m as 
measured by a nearby buoy. 



3. The PARFAIT approach 

In general, the altimetric information con- 
tent in the PARIS interferometric field phase 
will be very difficult to use. This is due to 
the impact of the incoherent component in the 
reflected signals. The incoherent component 
causes fading and winding. 

On the one hand, at a practical level, fad- 
ing events will prevent stable phase tracking of 
the complex field. Even if as in in the Bridge-2 
experiment the sea surface is relatively smooth 
and fading events are not so frequent, a single 
event can severly complicate the use of phase 
information if countermeasures are not taken. 
In general, however, the reflected field will fade 
very often. As discussed in |23j it is possible 
to inject in the system (during a fading event) 
a model-based phase to "glue" the phase his- 
tory, but this approach will in general necessi- 
tate the input of too much model information 
into the data in rough sea conditions. 

More importantly, as explained in the pre- 
vious section, the reflected field incoherent 
component will cause arbitrary winding of the 
field phasor. This means that the reflected 
unwrapped phase, unlike the direct one, can- 
not be directly used for ranging. Indeed, as we 
have shown in previous work I26j , the reflected 
field accumulated phase will generally wander 
around the complex plane, travelling to dif- 
ferent winding number kingdoms, even in the 
absence of fadings (see Figure EJ). That is to 
say, even if a very high SNR system is devised 
to get around the problem of field fadings, the 
interferometric unwrapped phase will not be 
directly usable for ranging. Unlike the prob- 

3 PARIS Interferometric Processor. 

4 PARFAIT stands for PARis Filtered-field Altimetric Tracking. 



lem of fadings, this is a fundamental issue, not 
a practical one. 

An approach discussed in |261 I19|. PIP 3 , 
involves the use of multiple frequencies for the 
synthesis of a longer wavelength which will be 
more immune to fadings. Here we discuss an- 
other approach, PARFAIT 4 , which is in fact 
complementary to PIP. 

In the PARFAIT approach, we begin by 
noting that although the reflected field un- 
wrapped phase carries no ranging informa- 
tion, this need not be a fundamental prob- 
lem. What is needed is the coherent geophys- 
ical field component in the signals, which is 
near zero frequency in comparison with the 
others — a sort of average field. This average 
field is just the coherent component in the 
reflected signals after downconversion. With 
this in mind, PARFAIT consists of the three 
steps described next. 

The first practical step to extract the co- 
herent part is to work with the interferometric 
field, the ratio for reflected versus direct com- 
plex field. This has the advantadge of error 
cancellation, e.g., in Doppler matching of the 
incoming signals, and of depending only on the 
lapse. 

The second step is to "counter-rotate" the 
interferometric field using an a-priori model of 
the reflection process. 

The third step is filtering the resulting 
counter-rotated interometric complex field to 
finally extract the coherent phase for esti- 
mation of the carrier lapse phase. Counter- 
rotation allows for longer filtering times. 
These are fundamental to extract the coherent 
component, which decays exponentially with 
the square of sea surface standard deviation 
(sea state) over effective wavelenght (the wave- 
length divided by the sine of satellite eleva- 
tion). 

Finally, the phase lapse information ob- 
tained from the couter-rotated, averaged, com- 
plex interfereometric field is used for altimetry. 

This new approach to PARIS altimetry is 
described in more detail in the following sec- 
tions. As we discuss, it has proven to be a 
very robust and precise processing method. 
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4. PARFAIT processing of the 
Bridge 2 dataset 

At low altitudes, simple geometrical con- 
siderations lead to the following equation re- 
lating the height of the receiver over the re- 
flecting surface (considering the same height 
for the upward looking antenna and for the 
downward looking antenna) with the lapse — 
the measured delay measured between the di- 
rect and reflected GNSS signals: 

C P {t) = cAt p (t) = 2h(t)-$m(e P (*)) + &, (6) 

where Cp(t) is the lapse in meters at time t, 
c is the speed of light, Arp(t) is the tempo- 
ral lapse in seconds, h(t) is the height of the 
bridge, ep(i) is the elevation of the GPS satel- 
lite with a specific PRN number P, and b is 
the hardware-induced delay bias, considered 
to be a common constant in time. A first esti- 
mation of the height of the receiver can easily 
be performed through a linear fit of the lapse 
with respect to the sine of the elevation angle 
of each satellite. 

In phase processing, the lapse is measured 
only up to an integer number of cycles N. 
Equation (j^J must be rewritten as follows 

Cp(t) = 2h(t) • sin (e p (t)) + b + N p X, (7) 

with Cp (t) is the carrier lapse in meters and A 
is the carrier (LI) wavelength. In other words, 
the equation for each satellite contains an ad- 
ditional unknown parameter, Np. In order to 
use all the satellites for one height estimation, 
it becomes necessary to estimate also Np, i.e. 
to solve the ambiguity problem. 

In order to solve the estimation problem, a 
minimization search is carried out for all these 
paremeters: h and b (as real constants) and 
Np (as integers). 

However, as discussed, the interferometric 
field should be first filtered to extract its co- 
herent component. Filtering should be long 
enough to extract the coherent component but 
short enough to keep the geophysical signals of 
interest pass through. 

This means that the geophysical coherent 
component we are after should not change for 
more than a small fraction of a cycle during 
the time duration of the filter. The maxi- 
mal allowable time thus depends on the ele- 
vation angle and rate of change of elevation 
of the satellite and, just slightly, on the tide 
motion. In the case of interest, it turns out 
the maximum filtering time should be around 



10 seconds. In other words, in 10 seconds, at 
least for one satellite, the coherent interfero- 
metric phase changes by more than 5 radians. 
With this filter length it is not possible to sep- 
arate coherent and incoherent components of 
the field, and fading events are not completely 
eliminated. However, a realistic estimation of 
the bridge height (and bias) does become pos- 
sible. 

As mentioned, to extract the coherent com- 
ponent a longer averaging period should be 
used. To this end, we first counter-rotate the 
interferometric field using a first guess of the 
bridge height, as we now explain in more de- 
tail. 

After downconversion and despreading, we 
can express the reflected complex field as a 
sum of the coherent and incoherent compo- 
nents, 

E{t) = A r e iC ^' x + 0{t). (8) 

Now consider that we have a first guess for the 
height and bias parameters, i.e., a model for 
the lapse C m . This model is used to counter- 
rotate the field: 

E cc (t) = E(t)/E m (t) (9) 

= A r e 2i5h (t) sin(e p (t))+i<5f> _|_ Q' 

Clearly, the phase of the coherent field in 
equation Q will now vary much slower than 
the phase of the original reflected field as a 
function of the elevation (and therefore time). 
This allows for a longer filtering time, and the 
exraction of the coherent component of the sig- 
nal (recall that 0{t) has zero mean). 

The equation which relates the countero- 
tated phase lapse between direct and reflected 
signal, the satellite elevations and the 5h (i.e. 
the error between the first guess of the bridge 
height and the real value) is 

Cp{t) = 2 6h(t) -sin (ep{t)) + 5N P X+5b. (10) 

This is the new equation to be used to fit the 
lapse versus sine of elevation straight line and 
infer the height offset of the bridge and bias 
(with respect to the first guess used to counter- 
rotate the field). 

In order to solve the ambiguity problem, 
a search is performed in the space of the in- 
teger n-tuples and the one that produces the 
linear fit with smallest residue is selected. It 
is important to point out that the n-tuple 
search space is drastically reduced by the prior 
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field counter-rotation. For example, if the 
guess is within ± half meter, the n-tuple sub- 
space to be scanned can be limited to those 
n-tuples whose components belong to the in- 
terval [—3, 3], centered on the first guess of the 
n-tuple, as obtained from a real (as opposed to 
integer) ambiguity resolution. 

Another way to reduce the cardinality of 
the subspace of the n-tuples to check is to con- 
sider that satellites with similar elevation an- 
gles cannot have very different integer ambi- 
guities. 

5. Results and comparison 

The PARFAIT algorithm described in the 
previous section, has been used to analyze the 
first 10 minutes of the Bridge-2 data, Part Al 
and to the first 10 minutes of Part A2. The fol- 
lowing steps have been performed accordingly 
in batches of 2 minutes: 

• The EM fields, direct and reflected, have 
been computed through the usual corre- 
lation process. 

• The complex interferometric field has 
been counter-rotated (Equation 0). 

• The counter-rotated field has been fil- 
tered using a 30 s window. 

• The phase of the interferometric, 
counter-rotated and filtered field has 
been unwrapped. 

• For every possible set of values of the 
ambiguities Np, a straight line has been 
interpolated to the phase histories (one 
for each visible satellite) against the el- 
evation angle (Equation (fTTTfl ^1 . The best 
fit has been identified. 

The analysis has been carried out for al- 
most 5 all visible satellites (see Table 2 and Fig- 
ure El)- The phase histories are shown in Fig- 
ure 11(a) A straight line has been fit through 



these phase histories, against the sine of the 
satellite elevation angle (Figure 11(b) I. 

This fitted line gives an estimation of the 
bridge height of 18.61 m, a hardware bias of 
-0.81 m and, as first guess for the n-tuple 
that solves the ambiguity problem, the val- 
ues [0 1 1 2 3]. Now, a search in a subset of 
T 6 is carried out to minimise the residuals of 
the fit in the space of the n-tuples of integers. 

Satellites outside the Zeeland Mask 1201 1211 are not considered (see also the caption in Figure H(H . 



The subspace considered is the one spanned 
by all the combination of integers between ±3 
around the first guess. The result is the n- 
tuple [0 2 2 4 5] which gives a bridge height 
estimation of 18.82 m and an instrumental bias 
of -0.45 m. 

This procedure has been implemented with 
data from the first 10 minutes of Part Al and 
A2 of the Bridge-2 experiment. The results 
are reported in Table 3 and in Figure ^| for 
Part Al and in Table 4 and in Figure ^] for 
Part A2. 

Fitting both parts to the tide curve, i.e. 
choosing the bias that minimizes the stan- 
dard deviation of the data to the available 
tide "ground truth" , leads to an altimetry 
bias of 40 cm and a standard deviation of less 
than 1 cm. This bias could have an origin in 
the ground "truth" , due either to an error in 
the determination of the absolute value of the 
height of the bridge performed using the GPS 
buoy available data (only a few seconds, which 
may have caused ambiguity resolution prob- 
lems) or, partially, to some anomalies in the 
water mass flow in the vicinity of the bridge 
structures. 

Moreover, considering also that the tide 
dynamics measured below the bridge may be 
delayed with respect to the place were the tide 
was measured, the best fit (over both bias and 
delay) is obtained with a delay of 1 minutes 
and 37 seconds with respect to the time of the 
tide data collection and with a bias of 40 cm. 
The standard deviation of the fitted data with 
respect to the tide curve is in this case of 0.3 
cm. 

To summarize, the proposed approach 
to PARIS altimetry, the PARFAIT technique, 
leads a very precise estimation of the tide, 

• without the need to insert any kind of 
model for the phase of the reflected sig- 
nal during fadings, 

• without rejecting too many visible satel- 
lites because of their poor SNR and/or 
frequent fadings. 

Finally, we note that this technique is di- 
rectly applicable for PARIS phase processing 
from air and spaceborne applications, as long 
as a suitable model for the lapse phase can be 
constructed. This will be the subject of future 
work. 
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PRN 


elevation 


mean S1MR 
(direct) 


mean SJMR 
(reflected) 


14 


17° 


29.4 dB 


25.0 dB 


25 


17° 


32.0 dB 


25.8 dB 


1 


30° 


31.2 dB 


24.6 dB 


7 


38° 


33.2 dB 


29.4 dB 


11 


62° 


34.0 dB 


29.4 dB 


20 


78° 


30.4 dB 


26.6 dB 



Tabic 2: Visible satellites, their elevation in degrees, the 10 ms coherent integration mean SNR^g, 
(20 log 1( ) [pcak-grass/grass correlation coefficient]) for the direct and the reflected signal. 



time 
(minutes 
from start) 


instrumental 
bias [cm] 


bridge height 
estimation [m] 


assumed 
ground truth 
[ml 


difference 
[cm] 


1 


-0.45 


18.83 


18.44 


39.71 


3 


-0.45 


18.82 


18.42 


39.63 


5 


-0.46 


18.81 


18.41 


40.19 


7 


-0.45 


18.79 


18.39 


40.21 


9 


-0.26 


18.78 


18.38 


39.83 



Table 3: Results of the bridge height estimation during the first 10 minutes of the Part Al data. 



Acknowledgements 

The authors wish to thank Manuel Martin- 
Neira (Technical Manager of the ESA/ESTEC 
Contract No. 14285/85/nl/pb under which 
this work was carried out) and Maria Bel- 
monte (ESA/ESTEC) for useful discussions 
and real collaboration. We also thank the 
other partners in the project, especially GMV 
for the GPS buoy data analysis. 

All Starlab authors have contributed signif- 
icantly; the Starlab author list has been ordered 
randomly. 

References 

[1] M. Armatys, A. Komjathy, P. Axelrad, 
and S. Katzberg. A comparion of GPS 
and scatterometr sensing of ocean wind 
speed and direction. In Proc. IEEE 
IGARSS, Honolulu, HA, 2000. 

[2] P. Beckmann and A. Spizzichino. The 

scattering of electromagnetic waves from 
rough surfaces. 1963. 

[3] M. Caparrini. Using reflected GNSS sig- 
nals to estimate surface features over wide 
ocean areas. Technical Report EWP 
2003, ESA report, December 1998. 



[4] M. Caparrini and G.Ruffini. Casablanca 
data processing. Starlab "Knowledge 
Nugget" kn-01 11-001, 2001. 

[5] E. Cardellach, J.M. Aparicio, A. Rius, 
J.S., and J. Torrobella. Application of 
the PARIS concept to transoceanic air- 
craft remote sensing. Technical report, 

2001. PARIS Alpha WP5200 - ESA Con- 
tract 14285/85/NL/PB. 

[6] E. Cardellach and A. Rius. Comparison 
of PAFEX estimates with ground truth. 
Technical report, IEEC, 2002. WP210 
of OPPSCAT 2 Project - ESA Contract 
RFQ/3-10120/01/NL/SF. 

[7] E. Cardellach and A. Rius. Inversion of 
PAFEX data with elfouhaily's technique. 
Technical report, IEEC, 2002. WP205 
of OPPSCAT 2 Project - ESA Contract 
RFQ/3-10120/01/NL/SF. 

[8] E. Cardellach and A. Rius. Preprocessing 
of PAFEX data. Technical report, IEEC, 

2002. WP120 of OPPSCAT 2 Project - 
ESA Contract RFQ/3-10120/01/NL/SF. 

[9] E. Cardellach, G. Rufflni, D. Pino, 
A. Rius, A. Komjathy, and J. Garrison. 
Mediterranean balloon experiment: GPS 
reflection for wind speed retrieval from 



2003 Workshop on Oceanography with GNSS Reflections 8 



time 
(minutes 
from start) 


instrumental 
bias [cm] 


bridge height 
estimation [m] 




assumed 
ground truth 
[ml 


difference 
[cm] 


1 


-0.27 


17.54 


17.11 


42.34 


3 


-0.28 


17.52 


17.08 


42.36 


5 


-0.26 


17.47 


17.04 


42.52 


7 


-0.08 


17.44 


17.02 


41.95 


9 


-0.08 


17.41 


16.98 


42.71 



Table 4: Results of the bridge height estimation during the first 10 minutes of the Part A2 data. 



the stratosphere. To appear in Remote 
Sensing of Environment, 2003. 

[10] J.L. Garrison, S. Katzberg, V. Zavorotny, 
and D. Masters. Comparison of sea sur- 
face wind speed estimates from reflected 
GPS signals with buoy measurements. 
In Proc. IEEE IGARSS, Honolulu, HA, 
2000. 

[11] J.L. Garrison, S.J. Katzberg, and C.T. 
Howell. Detection of ocean reflected GPS 
signals: theory and experiment. In IEEE 
Southeaston '97. IEEE, April 1997. 

[12] J.L. Garrison, G. Ruffini, A. Rius, 
E. Cardcllach, D. Masters, M. Armatys, 
and V.U. Zavorotny. Preliminary re- 
sults from the GPSR mediterranean bal- 
loon experiment (GPSR-MEBEX). In 
Proceedings of ERIM 2000, Charleston, 
South Carolina, USA, May 2000. 

[13] L. Garrison, S. Katzberg, and M. Hill. Ef- 
fect of sea roughness on bistatically scat- 
tered range coded signals from the GPS. 
Geophysical Research Letters, 25:2257- 
2260, 1998. 

[14] O. Germain and G. Ruffini. Least square 
inversion of PAFEX data. Technical re- 
port, Starlab, 2002. WP200 of OPP- 
SCAT 2 Project - ESA Contract RFQ/3- 
10120/01/NL/SF. 

[15] O. Germain, G. Ruffini, F. Soulat, M. Ca- 
parrini, B. Chapron, and P. Silvestrin. 
The GNSS-R Eddy Experiment II: L- 
band and optical speculometry for direc- 
tional sea-roughness retrieval from low al- 
titude aircraft. In Proceedings of the 2003 
Workshop on Oceanography with GNSS- 
R. Starlab, July 2003. 

[16] A. Komjathy. GPS surface reflection us- 
ing aircraft data: analysis and results. In 



Proceedings of the GPS surface reflection 
workshop. Goddard Space Flight Center, 
July 1998. 

[17] J. LaBrecque, S.T. Lowe, L.E. Young, 
E.R. Caro, L.J. Romans, and S.C. Wu. 
The first spaceborne observation of GPS 
signals reflected from the ocean surface. 
In Proceedings IDS workshop. JPL, De- 
cember 1998. 

[18] M. Martm-Neira, M. Caparrini, J. Font- 
Rossello, S. Lannclongue, and C. Serra. 
The PARIS concept: An experimental 
demonstration of sea surface altimetry us- 
ing GPS reflected signals. IEEE Transac- 
tions on Geoscience and Remote Sensing, 
39:142-150, 2001. 

[19] M. Martm-Neira, P. Colmenarejo, and 
G. Ruffini. Ocean altimetry interferomet- 
ric method and device using gnss signals, 
April 2003. U.S. Patent No. 6,559,165 B2. 

[20] M. Belmonte Rivas and M. Martm-Neira. 
GNSS reflections: First altimetry prod- 
ucts from bridge-2 field campaign. In Pro- 
ceedings of NAVITEc, 1st ESA Workshop 
on Satellite Navigation User Equipment 
Technology, pages 465-479. ESA, 2001. 

[21] M. Belmonte Rivas and M. Martm-Neira. 
GPS coherent reflections from a smooth 
marine surface. In Proceedings of the 2003 
Workshop on Oceanography with GNSS- 
R. ESA, July 2003. 

[22] G. Ruffini, M. Caparrini, and 
B. Chapron. Improved ocean and 
em models for in-silico spaceborne 
GNSS-R. Technical report, PARIS Beta 
WP3200 - ESA ESTEC Contract No. 
15083/01/NL/MM, 2001. 



2003 Workshop on Oceanography with GNSS Reflections 



9 



[23] 



G. Ruffini, M. Caparrini, O. Germain, 
F. Soulat, and J. Lutsko. Remote sens- 
ing of the ocean by bistatic radar observa- 
tions: a review. Technical report, PARIS 
Beta WP1000 - ESA ESTEC Contract 
No. 15083/01/NL/MM, 2001. Available 
online at http://starlab.es. 



[24] 



G. Ruffini, E.Cardcllach, A.Rius, and 
J.M.Aparicio. Remote sensing of 
the ocean by bistatic observations: 
a review. Technical report, OPP- 
SCAT WP1000 - ESA ESTEC Contract 
13461/99/NL/GD, 1999. Available online 
at http://starlab.es. 

Ruffini, O. Germain, F. Soulat, 
Taani, and M. Caparrini. GNSS-R: 
Operational applications. In Proceedings 
of the 2003 Workshop on Oceanography 
with GNSS-R. Starlab, July 2003. 

[26] G. Ruffini and 
intcrferomctric 

and experiment 



[25] G. 
M. 



F. Soulat. 
processor 



ysis 



PARIS 
anal- 
results, 



http://arxiv.org/physics/0011027. Tech- 
nical report, IEEC and GMV - ESA 
Contract 14071/99/NL/MM, 2000. 

[27] G. Ruffini, F. Soulat, M. Caparrini, and 
O. Germain. The GNSS-R Eddy Ex- 
periment I: altimetry from low altitude 
aircraft. In Proceedings of the 2003 
Workshop on Oceanography with GNSS- 
R. Starlab, July 2003. 

[28] V. Zavorotny and A. Voronovich. Scatter- 
ing of GPS signals from the ocean with 
wind remote sensing application. IEEE 
Trans. Geoscience and Remote Sensing, 
38(2):951-964, 2000. 

[29] C. Zuffada, R. Treuhaft, S. Lowe, G. Haij, 
M. Lough, L. Young, Wu S, and 
M. Smith. Altimetry with reflected GPS 
signals: results from a lakeside exper- 
iment. In Proceedings IGARSS 2000, 
2000. 



UL, - VIF 



Wtj + U)IF 



CO.s[(u) IF + bJd)t] 




) — ► 


LPF 






J L 





cos[((J d - u d )t] 



+ "IF 





LPF 




) ► 


J L 





sin[(LU IF + u>d)t] 



Figure 1: Starting from the in-phase component of the sampled signal, with a carrier frequency equal to the sum 
of the IF of the receiver and the Doppler frequency, a complex downconversion is performed. After low-pass 
filtering, the two obtained signals bring information about the amplitude of the backscattcrcd EM field for both 
the in-phase and in-quadrature components — the spread, complex field. 
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Figure 2: The histogram of the carrier phase variation, measured on an integration time interval (in this case 1 
ms). The x-axis represents cycles, while on the y-axis there are arbitrary units. The accumulation of 5-phase 
values around zero can be seen, as well as in the vicinity of ± half a cycle. 
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Figure 3: The directy signal carrier phase obtained accumulating the 5-phase according to Equation J^J. The 
stepped plot represents the accumulated phase as it is , i.e. without compensating for the navigation bit half- 
cycle variation, which is clearly visible. The lower curve represents the same phase after removal of this effect 
(navigation bit correction). The units on the x-axis are milliseconds, and on the y-axis thgy are cycles. 




Figure 4: Example of the tracked phase, without the Doppler contribution. The units are milliseconds on the 
x-axis and cycles on the y axis. The integration time is 20 ms. 
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(a) Direct field. 



(b) Reflected field. 



Figure 5: Example of field amplitude time series (phasor dust ball, Casablanca Experiment). The units are 
milliseconds on the x-axis and correlation coefficient units in dB on the y axis. The integration time is 20 ms. 
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(a) Direct field. 
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Figure 6: Example of complex field time series (phasor dusball, Casablanca Experiment). The units are 
correlation coefficient units on both axis. The coherent integration time is 20 ms. 
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Figure 7: Histogram of the phase for a noiseless simulated field. It can be observed that the unwrapped phase 
wanders around multiple winding number kingdoms, while tending to spend more time around an average 
complex field point. This illustrates the fact that the unwrapped phase cannot be used directly for altimetry, 
even in the absense of noise. 
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Figure 8: The simulated EM field at LI frequency (phasor dustball), after reflection on the sea surface. The 
simulation has been performed with the GRADAS software 1221 developed by Starlab. This simulation is for a 
wind speed of [/io = 3 m/s. 




Figure 9: The blue curve is the phase of the interferometric field for PRN number 7, from minute to minute 
8 of Part A2 of the Bridge 2 experiment. The occurrence of a fadings and of isolated cycle slips can be seen. 
These phenomena disappear in the phase of the filtered interferometric field (red line). 
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Figure 10: Each colored arc represents the position of a GPS satellite from the start of the Part Al of the 
experiment to the beginning of Part A2 plus 10 minutes. The green mask represent the area where the GPS 
signal reflections are supposed to be free of shadowing phenomena due to the bridge structure, and therefore 
only the satellite within this mask can be taken into consideration for PARIS processing. The bold parts of the 
lines represent the first and the second 10 minutes periods. 
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n (elevation angle) 



(a) In this figure, the reflected-minus-direct phase delay 
for each PRN is plotted with N P = 0. 
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(b) In this figure, the reflected-minus-direct phase delay 
for each PRN is plotted with N P = {001123}. 



Figure 11: Each colored spot represents the reflected-minus-direct phase delay versus satellite-elevation for a 
different satellite (PRN number). 
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minutes from start part At 

(a) The upper line is the bridge height over the water line 
estimated using PARFAIT. The bottom line is the bridge 
altitude according to the available tide measurements and 
the calibration of the absolute height using a GPS buoy 



minutes from start part At 



(b) In this plot, the estimates of the height (red *) and 
the actual height (line) are shown, after adding to the 
actual height the mean of the values of the last column 
of Table 3. 



Figure 12: The bridge height estimate during the first 10 minutes of the Part Al of the data. 
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(a) The upper line is the bridge height estimated using 
PARFAIT. The bottom line is the bridge height according 
to the available tide measurements and the GPS buoy 
measurement 1201 . 
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(b) In this plot, the estimates of the height (red *) and 
the actual height (line) are shown, after adding to the 
actual height the mean of the values of the last column 
of Table 4. 



Figure 13: The bridge height estimation during the first 10 minutes of the Part A2 of the data. 




Figure 14: The solid line is the distance between the up-looking antenna and the sea surface, according to the 
available tide measurements and the GPS buoy measurement 1201 . The green dots are the estimated values, 
after removing the vertical bias. 



